Two Anderson Impurities in the Kondo Limit. 
Systematic Study 



J. Simonin 

Centro Atomico Bariloche and Instituto Balseiro, 
8400 S.C. de Bariloche, Rio Negro, Argentina 
(Dated: today) 

We analyze the two Anderson impurity problem, in the strong Coulomb repulsion limit, by means 
of variational wave functions whose equations we solve analytically. We found two pairs of Doublet 
states, one odd and one even with respect to the midplane between the impurities. These Doublets 
make a significatively strong use of the hybridization terms of the Hamiltonian, and have a much 
larger Kondo- like correlation energy than a single impurity (the Kondo energy). Furthermore, 
these Doublets combine to form a Super-Singlet. This Super-Singlet makes use of the remaining 
hybridization probability to improve its energy with respect to the Doublets energy, but for low- 
Dimensional systems and at the inter-impurity distances where the Doublets are full developed, 
its energy gain is exponentially small. The interaction behind the Doublets is a first-order (in 
the effective Kondo coupling) one, it also generates a parallel alignment of the impurity-spins. This 
ferromagnetic impurity-impurity response is thus generated without resort to the RKKY interaction. 
The effective range of this first-order Doublet interaction is also much larger than the one of the 
RKKY interaction. 

PACS numbers: 73.23.-b, 72.15.Qm, 73.63.Kv, 72.10.Fk 



I. INTRODUCTION 

The two-impurity Anderson (TIA) model is a classic 
problem in condensed matter physics, where it repre- 
sents the interaction of two magnetic impurities embed- 
ded in a metallic hosli The Anderson Hamiltonian is 
a generic Hamiltonian that describes how localized or- 
bitals with a strong internal Coulomb repulsion behave 
when they are connected to/through a metallic bath. 
The current advance in nanotechnology allows for sev- 
eral man-made realizations of such systems^: Quantum 
Dots in the Coulomb blockade regimen connected via a 
low-dimensional electron gas^, magnetic atoms individu- 
ally deposited on metallic surfaces, nanoclusters or mag- 
netic impurities connected to/through metallic carbon- 
nanotubes, etc. This kind of circuits, characterized by 
the low-dimension of the metallic host, promises to be a 
relevant circuit component in quantum electronics. 

The TIA Hamiltonian has been the subject of many 
theoretical studies, ranging from perturbation theoryi 5 ^ 
and narrow band approximation 7 to renormalization 
group analysis&SiiSiiiiiSii&ii, conformal field theory and 
Quantum Monte Carlo simulation*^. With some 
approximations/simplifications, as detailed in Ref.pi], 
some analytical results has been obtained. An other the- 
oretical approach to Anderson Impurity systems has been 
the use of variational wave functions (VWF)A 7 ^i£*i2i2iI. In 
Ref. a set of VWF equations that describes the lower 
energy state of the TIA system has been numerically an- 
alyzed. In this paper we present a systematic VWF anal- 
ysis of the TIA in the strong-Coulomb-repulsion Kondo 
Limit. By means of VWF we obtain several set of equa- 
tions that describe the lower energy states of different 
subspaces of the system. We analytically solve those 
equations, identifying two Kondo-like energy scales 9,13 



that depend strongly in the inter-impurity distance. The 
higher one corresponds to the formation of Kondo-like 
Doublet states and the lower one corresponds to the for- 
mation of a "composite" Kondo Singlet, based upon the 
screening of the Doublet states. Our equations depend 
on the details of the metallic host through a simple Co- 
herence factor that, when needed, we evaluate in the one 
Dimensional (ID) case because of its importance for tech- 
nological applications. The 2D and 3D scenarios are also 
analyzed. 

We found that the Doublets states generates a ferro- 
magnetic impurity spin-spin correlation without resort to 
the RKKY interaction. The effect of the RKKY terms 
in the Doublets structure and energy is also analyzed. 

Whereas our analytical results confirm the general pic- 
ture of the system behavior that is already known from 
numerical simulations, they also point out that the main 
interaction between the impurities, contrary to the gen- 
eral beliefSi*2S, is not given by the RKKY process but by 
the simplest electronic process behind the formation of 
the Kondo-Doublet states. 

This paper is organized as follows: In Section [n] we 
write down the TIA Hamiltonian and we make a basis 
change to the mirror symmetric basis. In Section IlIII we 
make a preliminary study of the system, analyzing the 
main features of the Hamiltonian at some special dis- 
tances between the impurities. In Section IIVI we make 
a systematic analysis of the different subspaces of the 
Hamiltonian. The effects of the intrinsic RKKY inter- 
action in our states are analyzed in Section TlV Dl The 
high-Dimension scenarios are analyzed in Section[V] The 
Impurity spin-spin correlation is evaluated in Sect ion Ivfl 
In Section IVIII we discuss the conclusions and perspec- 
tives open by the method, its possible extension to simi- 
lar situations and its relation to previous work. 
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II. THE HAMILTONIAN 

The Anderson Hamiltonian for magnetic impurities di- 
luted in a metallic host is: 



ko 



jka 



-E d J24 a d j(T + UJ24A d ^d j l , (1) 

j 

where c^a and dj a are fermion operators which act re- 
spectively on the conduction band states and on the or- 
bital states of the magnetic impurity placed at rj . Single 
state energies e k ,—Ed are referred to the Fermi energy 
(Ep = k F /2m), i.e. there is an implicit — pL N term in 
the Hamiltonian that regulates the population of the sys- 
tem {fi = Ep is the chemical potential and N the total 
number operator), V is the d— c hybridization and N c is 
the number of cells in the metal. In the Kondo limit that 
we analyze here the impurity levels are well below the 
Fermi energy (— Ed <C 0), and can not be doubly occu- 
pied due to the Coulomb repulsion in them (U — > oo). In 
order to simplify calculations we renormalize the vacuum 
(denoted by \F)) to be the conduction band filled up to 
the Fermi energy and we make an electron-hole transfor- 
mation for band states below the Fermi level: bl = c k a 
for |fc| < kp. In this way the energy of a hole excitation 
is explicitly positive. 

We consider here two impurities placed one at — r/2 
(the Left impurity) and the other at r/2 (the Right im- 
purity). 

In this situation (—Ed < 0, U — > oo) the lower en- 
ergy configurations of the system are the ferromagnetic 
(ferro, FM) triplet ( d^d^ \F), ... ) and the antiferro- 

magnetic (antiferro, AF) singlet ( [d\ j ^d* R ,— rl ' 
with energies —2Ed- These states are not eigenstates of 
the system because the hybridization terms (Hy) mixes 
these states with configurations having excitations in the 
band. 

An important state for magnetic impurities in a metal 
is the Kondo singlet, the ground state when a single im- 
purity is considered. If this impurity is at the origin, this 
singlet is described by the VWF— 

\ s ») = \ F )-J2 z ® ( b \Ai\ F ) +b \i d U F )) > ( 2 ) 



where the function Z(q) is to be determined variationally. 
It turns out that Z(q) — v/{E K + E d — e q ). A self- 
consistent equation is obtained for the singlet energy Ek, 
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Ek + Ed — e q 



(3) 



where the q sum is over hole excitations (\q\ < /cp), and 
v = V/^/W c . After the form 



is assumed Eq.Q becomes an equation for the Kondo 
energy gain S K , 



£ d + ( 5 A . = 2v 2 ]r— 1 



5k + e Q 



(5) 



From here, in the E d , D 3> 8k regime, it is obtained 



5k = D exp - 



1 

2 J n. 



(6) 



In this equation J„ = n V 2 / Ed , n being the density of 
states at the Fermi level and D the half band- width. J n , 
the effective Kondo coupling, is the relevant parameter 
for these systems, and D gives the energy scale. If we 
demands that 5k be one thousandth of D, a reasonable 
value, it follows J n = 0.07238. 

Note that although all the b^ qa .d L -\F) configurations 
enter in the definition of the Kondo singlet only the holes 
with e q < 5k make a significant contribution given that 
the Z(q) variational amplitude results proportional to 
1/(5k + e g ). The Fermi state \F) plays the role of a 
nearly virtual state in the Kondo singlet. 

The Hamiltonian Eq. conserves the total number of 
electrons Nt, the total spin St, its projection S z , and 
has also a mirror symmetry plane at the middle point 
between the impurities. For example, the FM state has 
Nt = 2 (plus a half filled conduction band), St,S z — 
1, 1 and is odd under the mirror reflection (d^d^lF) = 

- d R1 d Li\ F )), the AF state has Nt — 2 , St — and 
is even under reflection. Assuming a canonical ensemble 
in the thermodynamic limit, this allows the Hamiltonian 
to be analyzed by subspaces corresponding to different 
quantum numbers. 

It is usefully to rewrite the Hamiltonian Eq. in terms 
of one-electron operators of well defined mirror symme- 
try. This is easily done by changing to the symmetric 
and antisymmetric combinations 

1 , X 1 / \ 

CSka = -7^{C ka + CfcJ, CAka = ~J={Cka ~ CjfJ, 



V2' 



V2' 



dsa = ~^^ Ra + d La), dAa = ^(^^ _ dha) , (7) 

where k stands for (—k x ,k v ,k z ). The impurities are on 
the x-axis, at x = ±r/2 and the mirror plane is the x = 
plane. The transformed Hamiltonian is 

H = E e k J xka c xk a -E d ^ d\ a d X a + 



Xka 



Xc 



E K 



-E d 



(4) 



2v 5Z( C0S ~T" d Sa c Ska +isin^- d\ a C A ka + h.c), (8) 



plus eight four-<i-operator U Coulomb terms, which can 
be evaluated when needed. In the Hamiltonian above 
X = {A, S} and the fc-sum runs only in the right k- 
halfspace (k x > 0). The electron- hole transformation for 
the conduction band states in this basis is just b^ xka = 
c X kw for |fc| < kp. 
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The effect of the U terms in the U — > oo limit is to 
inhibit two of the six two-electrons in the two-impurity 
states. The forbidden configurations are 

(4|4t +4;4t)D and ( d Al d U ~4t4j,)l) ' 

i.e. the symmetric and antisymmetric combinations of 
the forbidden configurations in the direct basis, ^j^Lf!) 
and 4i4tD • The four allowed configurations are 

|AF) = (4 ; 4 T - 4^^)!), |FMt) - 4 t 4tI>- 

ifmo) = (4i4r +4A)i>. W) = < MJ>» ( 10 ) 

the first one is the AF two-impurity configuration, an 
even state, the last three are the S z = 1, and —1 com- 
ponents of the FM triplet. Higher impurity population 
states are forbidden. 



III. THE MODEL 

We base our theory on a careful analysis of the Hamil- 
tonian in the symmetrized basis, Eq.JjJJ. First we notice 
that in this basis the hybridization terms are decoupled, 
i.e. symmetric ds impurity orbitals only mix with sym- 
metric cs conduction band states and the same is true 
for the antisymmetric ones. Thus it is appropriate to 
think in terms of two nearly independent magnetic impu- 
rities, the symmetric impurity (SI) and the antisymmet- 
ric impurity (AI), keeping in mind that these are linear 
combinations of the two original impurities. The sym- 
metrized impurities are correlated through the action of 
the Coulomb interaction because the forbidden impurity 
occupations are not the double occupied SI (or AI) config- 
urations ( 4l4tD and 4l4ttl) ) the combinations 
thereof given in Eq.|0. 

Secondly, we notice the factor of 2 in the hybridization 
term of Eq.©, which originates in the change in the 
matrix elements of the Hamiltonian due to the change of 
basis states. This has important consequences because 
the correlation energy for a magnetic impurity Eqs.JS] 
EJ) depends on the effective hybridization raised to the 
power two. 

A. A singlet 

As a consequence of the above, when the two impurities 
are at a distance such that | cos ~ 1 and given 

that only the conduction holes near the Fermi level are 
relevant for the formation of a Kondo singlet, the effec- 
tive coupling for the SI impurity is nearly twice the bare 
coupling V whereas the effective coupling for the AI im- 
purity is near zero. Therefore, the maximum Kondo-likc 
correlation energy 84 that can be obtained by forming a 
SI Kondo singlet is given by 

5 4 = D exp- — . (11) 



Notice the factor 4 in the exponential, while the single 
impurity Kondo energy 8 k has only a factor of 2. The 
effective hybridization coupling provides a 2 2 factor but 
half the g-states in the sum of Eq. JSJ) have already been 
summed in doing the basis change, and this is the reason 
for the change from 2 to 4 in going from 8 k to 54. This 
result is exact for r = and a good approximation for 
the next few maxima of the effective SI hybridization 
coupling (r ~ n\p, n a small integer, the Fermi wave 
length) . This is a very large correlation energy gain for 
the two impurity case compared with the single impurity 
Kondo energy. For the above given value of J n — 0.07238 
, and assuming D = Ep = 10, 000K, one obtains 8k — 
10K and 84 = 316-KT, an remarkable effect. In fact the 
84/ D ratio is given by the square root of the 8k /D ratio, 
thus this magnification effect is higher the lower 8k /D. 

Let us examine this situation in more detail. Consider 
first the values of r for which the effective symmetric hy- 
bridization coupling is enhanced and the antisymmetric 
one is depressed, | cos (kpr/2)\ ~ 1, | sin (kpr/2)\ ~ 0, i.e. 
r ~ n\p. To analyze these regions we can write down 
the two-impurity SI singlet \SI), changing b\ cr d\_ C! by 

"Sqo d s—<T m th° Kondo singlet Eq.Q. The maximum 
correlation energy gain for such an state is £4, thus the 
energy of the SI singlet is Esi > — Ed — 84 , i.e. a very 
high correlation energy gain but a not too low total en- 
ergy because the average total population of the impuri- 
ties is 1 for such an state. In the appendix ^ we show 
that Esi is always greater than the ground state energy 
expected when the two impurities are far apart, 

E 2K = 2E K = -2E d - 28 K , (12) 

the ground state energy for the two magnetic impurity 
problem in the r — > 00 uncorrelated limit. 



B. A doublet 

We can try to take advantage of the correlation en- 
ergy gain 84 of the \SI) singlet by making the SI singlet 
the core of an odd doublet formed by adding an elec- 
tron in the inactive AI impurity, i.e. forming the doublet 
I -Dot) = d At\SI)- This state has an average total popu- 
lation of 2 electrons in the impurities but it contains the 
4t4il) impurity configuration, which is not within the 
configurations allowed by the U terms of the Hamilto- 
nian, Eas. (|9ll0|l . This can be corrected by noting that 

2 4 T 4j) = (4 T 4i + 4i4 T )i) + 

(4 T 4i-4i4 T )i) ■ ( 13 ) 

The first combination is the |FM0) allowed state and the 
second a forbidden one, thus we project it out of the 
wave function. In doing this we are loosing 1/4 of the 
active states of the \SI) singlet, we will see that this lose 
reduce the factor 4 of Eq. (|11|) to a maximum of 3 for the 
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correlation energy of the doublet. Thus, the odd doublet 
is given by 



|A, T )=c4 T |F) + 



S k {2Z 1 (k) &LJFMT) 



z 2 (k) &L T i FM °)} ( u ) 



where Sk = cos ^g 1 . We use the notation of Eq. (jTU|) 
for the doubly occupied impurity states. We apply the 
standard analytical Euler/Lagrange minimization pro- 
cedure to the expectation value of the energy , Ex — 
(X\H\X) / (X\X) , of this state. For the variational func- 
tions in the amplitudes of the wave function components 
we obtain 



Z 1 (k) = Z 2 (k) 



Er 



2E d 



(15) 



We see that they depend on the energy of the doublet 
Ed o , which is determined by the selfconsistent relation 



Er>„ — —Ed 



1 + cos q x r 
E Do + 2E d -< 



(16) 



where the sum index runs over the symmetric hole states, 
i.e. \q\ < kp and q x > 0. As we are looking for the 
lowest energy solution of this equation, we assume Eu = 
—2E d — 8 - Ea. i|16|) becomes a selfconsistent equation for 

So, 



3v 



2 1 + COS Qx T 

Y *o + e g 



(17) 



Two limits can be immediately obtained by just compar- 
ing this equation with the selfconsistent equation for Sk , 
Eq.JSJ. At r = one has cosg^r = 1 for all q's, thus 

5 3 = S o (r = 0) = D exp--i- . (18) 

o J n 

This is still a great enhancement. For the values used 
in the previous analysis we obtain £3 = 100K (whereas 
Sk = 1QK). This is an exact result because the AI im- 
purity is strictly inactive for r — 0, this result has been 
previously found in Ref.fS^]. For r>0 the contribution 
from the cos q x r term in the sum vanishes by decoherence 
effects, therefore S a (r > 0) = D cxp(-l/(| J n )) < 5 K - 
This is a physically wrong result, one expects S = Sk 
in this limit, but by the same decoherence effects the AI 
impurity is active in this limit even if r — MXf , M being 
a very big integer, and its contribution to the odd dou- 
blet can no longer be ignored. We correct this point in 
Section |IVU 

For the values of r other than the analyzed limits, we 
must evaluate the sums in Eci. (|17|) . The first one is the 
Kondo integral 



Ik{8) = 



1 



n Q N c ^ S 

q 



= ln( 



D 



(19) 



and the second one can also be summed analytically. We 
call it the Quantum Coherence Integral, and in ID we 
obtain 



n N, 
,f>S, 



E 



cos q x r 
S + e a 



(20) 



cos( PX )[Ci(px)-Ci(^)] 



rPS, 



sm(px)[Si(px) - Si ("7j)] . 



where p = kpr, \ =^rT an d Ci (Si) is the Coslnte- 
gral (Sinlntegral) function, as defined in Mathematica®. 
The Quantum Coherence Integral has a logarithmic de- 
pendence on S that goes like Ik ■ Notice that Iq (S, 0) = 
Ik(S), so it is useful to define the two- impurity Coherence 
factor Cq{5,t) = Iq(8,t)/Ik{&), which is a decaying os- 
cillatory function that depends weakly on S. In FigE]wc 
plot Cq as a function of r for various values of 8, the in- 
dependence of the 1D-RKKY interaction, Jr{t) : is also 
plotted for comparison. It can be seen that Jr{t) de- 
cays rather rapidly, with range r ~ Xf, whereas Cq has 
the much larger range ( = Xf Ef/Sk ), the Kondo 
length. In ID the dominant decoherence effect is that of 
the energy width ( S) of the packet of holes that forms 
the Kondo cloud. In higher Dimensions they appear also 
angular effects that contributes to the decoherence with 
typical length of the order of Af- 
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FIG. 1: The two-impurity ID Coherence factor, Cq(S, r) and 
the r-dependence Jr(t) of the RKKY interaction as func- 
tion of r. See that Cq depends weakly on S. The period of 
Jn(r) is half the one of Cq and its amplitude decays rather 
rapidly. These characteristics are due to the fact that the 
RKKY depends on two excitations, whereas the Coherence 
factor involves just one. 

Using the results of the sums Eqs.(Hni EDI, Ea. l(T7|) 

reads 

E d + 5 = V 2 n l(l + C Q (S ,r))\n^r^, (21) 
z do 



which can be recast as 



S a (r) = D exp - 



\[l + C Q (S ,r)]J n 



(22) 
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still a selfconsistent equation but given the weak depen- 
dence of Cq on S, it converges very rapidly when used 
iteratively. It can be easily checked that it gives the 
two limits already discussed, corresponding to Cq = 1 
(r — > 0) and Cq — (r — + oo). We will show later, by 
evaluating the energy of the full odd doublet, that this 
approximation is in fact good enough at r ~ u\f, more 
precisely, around the maxima of Cq. At r ~ (n 4- 1/2) Af 
the more active impurity is the AI one. Therefore for 
these points we could repeat all the arguments of this 
section but for an even doublet based on the AI singlet, 
however it is enough to observe that the square of the 
effective hybridization coupling of the AI impurity goes 
as sin(fc x r/2) oc (1 — cos (k x r)). The final result for 
the correlation energy gain of the even doublet is in fact 
Eq. (|22|l but with a minus sign in front of Cq . We show 
these energies in Fig|5]as a function of r. We see that be- 
tween two consecutive maxima of each of them there is a 
flat region in which the other one has a maximum. This 
coincides with the behavior of (1 ± Cq(t)) exponentially 
amplified. See also that when S is at one maximum S e is 
very small and viceversa, for the first maxima their ratio 
is exponentially small. 
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ground state of the SI impurity. Furthermore, this dou- 
blet, at r = 0, is probably the ground state of the system, 
given that coherence effects are at their maximum and 
the AI hybridization is strictly zero. Furthermore, when 
5 a is at one of its maxima, there is very little energy to be 
gain from the AI impurity, as we can see from de value of 
5 e at those same points. Thus it seems that it make little 
sense to modify further the structure of these doublets. 

However, we see that the vertex state of the \D a ^) dou- 
blet is the d A ^\F) configuration, and d A , \F) the one cor- 
responding to the S z = —1/2 component of the doublet. 
These states are the ones that can be connected through 
an AI singlet. Therefore we try the following supersinglet 

\SS) = \F)-zJ2 Z{k) A k {b\ ki d% + b\ k ^ Ai )\F) 

k 

-iJ2Y(k,q) A k S q { {b Aki bl q] +b\ k] bl qi )\YMQ) 

k . q 

+ 2 6 L/s<J FM T> + 2 b Ak / Sql \FMl)} , (23) 

where A k — sin^J^. The variational amplitudes Z(k), 
Y(k, q) and the supersinglet energy Eg turn out to be 

Y{k, q) = v Z(k)/(E S + 2E d -e k - e q ) , (24) 
Z(k) = 2v/D z (k), (25) 
E s = 2 v^(l-cos(fe x r))Z(fe) , (26) 



where the denominator of Z{k) is 



D z {k) =E s +E d -e k -3v 2 



1 + cos q x r 



^ Eg + 2Ed — e k 



To find Es we must transform this denominator. It 
closely resembles Ea. (11611 for Ed o , with the change 
Eb c i—* Es — e k . In fact D z (k) comes from the contribu- 
tion of the odd doublet to the \SS) energy. Therefore we 



FIG. 2: Correlation energy gains of the odd 8 and even S e 
doublets as a function of the distance between the impurities, 
for J n = 0.07238 (S K = 0.00LD). 

From Eas. l|17l22|l . and its equivalent for the even dou- 
blet, one can say that (1 + Cq{8d 1 t))/2 determines the 
relative strength of the square of the SI impurity hy- 
bridization as function of r better than cos (k F r/2y. In 
fact Cq (5d , r) has embodied in it the r decoherent ef- 
fects on the packet of the u 5d relevant holes of this 
two-impurity Kondo-like interaction. 



C. A supersinglet 

At this point of our analysis we see that the odd dou- 
blet seems to obtain all the Kondo like correlation energy 
of the SI impurity. This is because it is formed out of the 



Ec 



E 



Do -7 



2E d - 6 



(27) 



Thus, using {Ik, Iq, Cq } for the sums in D z , and then 
Ea.l(?I)l for E d + S a , -D z (k) becomes 



D z (k) 



-V z n {l + C Q {5 ,r))\n 



(7 + e/c) + 

So + 7 + e k 



(28) 



now, note that if 7 > 0, i.e. Es < Ejj o , —D z (k) is always 
positive. Furthermore, we expect this reduced version of 
the supersinglet to be good just in the regions in which 
the odd doublet is. In those regions the contribution 
to be expected from the AI hybridization channel is very 
small, i.e. S 3> 7, and thus the holes that will contribute 
the most to the sum in the Eg equation Ea. (|26fl are the 
ones at an energy distance of the order of 7 from de Fermi 
level. In these conditions the argument of the logarithm 
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in Ea. (|28(l is very close to one, and thus is safe to take 
—D z (k) = (7 + e i- ) when used in the sum of Eg . (|26[) . The 
next order correction (ln(l + x) ~ x) can be included 
at no mathematical cost, but it requires introducing a 
new parameter, the Ed/ D ratio. We keep it simple just 
disregarding this correction. Therefore, Eq. (|26fl becomes 



2 E d + 8 + 7 = 4v 2 



(1 - cos (k x r)) 



1 + e-k 



2V 2 n [l-C Q ( 7) r)} In 



from which 



7(7") = D exp - 



[l-C Q ( 7 ,r)]J n 



(29) 



(30) 



This is a very small number when S (r) is strong (Cq 
near 1). Notice that the Cq factor in Ea. (|30[l depends 
on 7, not on 8 , thus it is closer to 1 than Cq(<5 ). 

Near the second maximum of 8 Q , at r = I.IAf, and 
for J„ = 0.07238 Ea. i|22[l converges in five iterations to 
S = 0.00324 (and Cq = 0.607) starting with 8 Q = 8 K as 
iteration seed. Ea. (|30|l converges in seven iterations to 
7 = 1.394 10~ 29 (Cq = 0.792). Near the first maximum, 
at r = 0.1A F we obtain 8 a = 0.00922, 7 = 5.70 10~ 124 . 

These are very small numbers indeed, but 7 is positive 
and small enough to justify the approximations made 
to reduce Eq.J^HJ, and thus the supersinglet is a better 
option as a possible ground state of the system than the 
odd doublet. At r = 0, Cq = 1 and 7 = 0, thus E s = 
Ed and the supersinglet collapses into the odd doublet, 



it is an uncorrelated combination of the b Ak 



and the b 



\D„ 



J Ak F [\ ±J °M states. At r — > 00 (Cq = 0), or at 
the Cq(t) = points, this version of the supersinglet 
fails to give some expected results, mainly that it must 
be E$ — —2Ed — 28 k for r — > 00. This is not surprising 
since we already expect this "partial" supersinglet not to 
give good results outside the r ~ n\p regions. 

The previous analysis give us a clear indication of what 
to look for in the different subspaces of the Hamiltonian, 
the chemical potential tell us no to go too far from Nt — 
(plus the half filled conduction band) thus, given the 
limits that U impose in the population of the impurities, 
Nt = 2 seems to be a reasonable upper value. In the 
next section we proceed to a systematic analysis of those 
subspaces, looking for the lowest energy state of each 
one of them. We start with the Nt — 1 one and finish 
with the one corresponding to Nt = 2. This order is 
not arbitrary because the Nt — 1 subspace is the less 
explored of them and Nt — 2 the most explored one. 



IV. SYSTEMATIC ANALYSIS 

A. Doublets subspace, Nt = 1 

The Nt — 1, odd subspace is the odd doublet sub- 
space. As follows from the previous analysis the odd dou- 



blet lacks the action of the AI hybridization channel on its 
vertex state, i.e the b\ 



>d\ t d\AF) and b\, r ,d\ 



jt 



J Akt u, Al u A1\- L ' I "Akl^A^A] \F) 

configuration groups. The last group is forbidden by 
Pauli exclusion principle and must thus be ignored. For 
the first group we must use the same procedure as in 
Eq. (|13|) . Thus we analyze the following variational wave 
function in this subspace 

I £> oT ) = dl r \F) + 1 Y, Z A (k)A k b\ k] \AY) 

k 

+ Y / Zs(k)S k {2bl kl \FM^+b^\FM0)} , (31) 

k 

the now full odd doublet. The selfconsistent equations 
for the variational amplitudes and energy are 



Z A (k) = Z s {k) = 
= -E d + 2v 2 Y 



Ed + 2Ed — efe ' 



E Do + 2E d 



(32) 
(33) 



After expanding the S 2 , A 2 , factors in Eq. H33|) . and as- 
suming Etj„ = —2E d — 5 a , Eq. (|33|l transforms into 



2 + cos q x r 
S + e„ 



(34) 



now, using { I K (S), and C Q (S,r), Eas. fiMty }, Eq.® 

becomes 



S a (r) = D exp 



-1 



[2 + C Q (S ,r)} J n ' 



(35) 



in the D,E d 8 limit. This formula for 8 (r) re- 
places the previous partial result Ea. (|22|l . It can be 
taken as an approximate solution by evaluating Cq at 
8k, or used iteratively to find the exact result with fast 
convergency. At its first maxima, i.e. at { Cq ~ 1, 
r ~ n\p }, 8 a < 83 as previously analyzed. At the 
Cq{8, r) = points, including the r — » 00 limit, Eq.^ 
gives S o {r / C Q (r ) = 0) = 8 K (= 8 2 ), the physically 
expected result. 

In the Nt = 1, even subspace, the relevant state is 
the even doublet, its construction process is the same 
as for the odd doublet but taking dg^\F) as the vertex 
state. The definition of the even doublet is similar to 
that of the odd doublet, Eq. l|31|l . but with the change 
A <-> S in the creation operators and in the hybridization 
amplitude factors A k ,Sk- A global minus sign appears 
for the Zx(k), given that for the AF, FM states we use 
the d) operators ordering shown in Ea. l|10|) . the one of 
the allowed impurity configurations. The energy Eo e of 
the even doublet states is given by 



E Dc {r) = -2E d -S e (r) 
8 e (r) — D exp ■ 



[2-C Q (8 e ,r)} J n 



(36) 



If we compare Ea. (|36|l for 8 e (r) with Ea. l|35|) for 8 Q (r) , we 
see that the sign in front of Cq changes, this is because 
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of the A k <-> S k change, i.e. the SI and AI impurities 
interchange roles. The analysis of 5 e as a function of r is 
also very similar to that of 5 (r) . The odd doublet corre- 
lation energy gain S (r) is higher than the even doublet 
gain 5 e when Cq{5k,t) > 0, and the opposite is true for 
Cq{5k,t) < 0. They are equal, and equal to Sk — &2, at 
the tk distances for which Cq(5k, Tk) = 0. 

These correlation energy gains of the odd and even 
doublets, S and 5 e , are plotted in Fig as a function of 
r for Sk — 0.001 D. For very large r, such that Cq(Sk, r) 
vanishes, the energy gain of the doublets tends to that 
of one Kondo singlet. This is not surprising since the 
doublets, in this limit, can be written as 

\Dxt) = \Li)®\S R )±\S L )®\Ri), (37) 

i.e., the combination of a Kondo singlet in one of the orig- 
inal impurities and the other single occupied plus (minus) 
the mirror image of that state. The total energy of this 
state is —2E,i — 8k, as for the doublets in the uncorrelated 
limit (r — > oo). 




FIG. 3: ID Correlation energy gains of the odd S and even 8 e 
doublets as a function of the distance between the impurities, 
for J n = 0.07238 {8 K = 0.00LD). The inset is 5 e up to very 
long distances, it periodically beats 25k up to r ~ £k- 

These energies gains must be compared with twice the 
Kondo energy, because, for very large r, the latter is 
the energy gain for a state that has simultaneously, and 
in a uncorrelated way, both impurities forming a Kondo 
singlet (\Sl) <8 \Sr)). The doublet gains So alternatively 
surpass 25k up to r ~ 

This effect, <5> = m&x(8 ,8 e ) > 25k, is present for 
values of J„ lower than J„ c = 0.240 (defined by 63 = 
2^2); the lower J n , the greater this effect. The periodicity 
with which \D a ), \D e ) alternate as ground states of the 
Nt = 1 subspace is twice the RKKY period (they are 
ruled by Cq{5k)), and they stand in this situation up to 
very long distances. These properties are due to the fact 
that the doublet Coherence Integral involves just one hole 
excitation. At the intersection points one has 5 e ,5 = 5k 
and Cq = . 



The electronic process behind these correlated Dou- 
blets is second order in the hybridization coupling V. 
The jump of a band electron from below the Fermi level 
to the impurities and then back to its initial place is the 
only process allowed by the Doublet VWF Eq. [pIT|) . 

B. SuperSinglet, 7V T = 

In the Nt = subspace there is the supersinglet, its 
full expression, taking into account both the even and 
odd doublets components, is 

\SS) = \F) - £{[ * Z (k) A k {b\ ki d% + b\ k ^ Ai ) 

k 

+ Z E (k) s k (&t 

^ST + ^Sk^Sl) 

k , q 

{Y (k,q)A k A q b\ ki b\ q] + Y E (k,q)S k S q tf Ski bl qr ) |AF) 
+ iY(k, q )S k A q [{b\ qi 

h \k\ + b Aq1 b Skl) l FM 0) 

+ 2 b AAi i™t> + 2 &Vk T |FMi) ]} • (38) 

The Zo(E)(k) configurations are the odd (even) doublet 
component of the supersinglet. The "origin" of each one 
of its components can also be traced down looking at the 
A k , S k factors and/or the results obtained for their vari- 
ational amplitude functions. After the functional mini- 
mization procedure, these functions turn out to be 

Y (k,q) = v (Z (k) + Z (q)) /D Y (k,q) , 
Y E (k,q) = v (Z E (k) + Z E (q)) /D Y (k,q) , 
Y(k, q) = v (Z (k) + Z E {q)) I Z?y (fc, q) , (39) 

where Dy{k, q) — {Eg + 2Ed — e k — e q ) and 

Z Q (k) = -[ 2v + v 2 - cqr) Z (q)/D Y (k, q) 

q 

+3v 2 ^(l + ccp)Z E {q)/D Y (k, q) } I D Zo (k) , (40) 
1 

D Zo (k) = -E S -E d + e k + 2v 2 £ , (41) 

where cqr = cos q x r. For Z E {k) a symmetric expression 
holds, with the changes, in Eas. l|40l41|l . E «-* 0,A <-> 
S, cqr 1 ► —cqr. For the supersinglet energy Es, the 
following expression is obtained 

E S = 2vJ2l(l-ckr)Z (k) + (l + ckr)Z E (k)} . (42) 

k 

To solve these equations, which are self consistent and 
recursive, one must proceed in a similar way as in the 
case of the partial supersinglet analyzed in Section UlI CI 
We work them out in the same approximation as before 
disregarding higher order terms that depend on the E^/D 
ratio. 
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First, we analyze the points at which coherence effects 
are suppressed, i.e. r — * oo and the tk points such that 
Cq(5k,tk) — 0. For these points 5 Q = 5 e = 5k and 
thus, assuming E$ = —2Ed — 5k — 7, we obtain 



Z x = -2v/D z 
Dz x = (7 + efe) 

2E d + 5 K +l = 8v 2 ^ 



1 



7 + efc 



(43) 
(44) 

(45) 



from which 7 = 5k follows. Therefore, the energy of the 
supersinglet at these points is 



Ec 



-2E r] - 25 



K 



(46) 



The value Cq{5k, tk) — marks the boundary between 
two regions. For Cq{5k,t) > the odd doublet is the 
lower energy doublet and it is the main component of the 
supersinglet, conversely, for Cq(5k,t) < the behavior 
of the supersinglet is determined by the even doublet. 
We work out now the energy of the supersinglet in the 
odd doublet region. In these regions 5 > 5 e and thus 
we take E$ = —2Ed — 5 — 7. With this replacement 
Eas. l|4UI42|l become 



Dzo 
E s 



(7 + e*), D Ze = (-f + A s + e k ) 
ckr 1 + ckr , 



D Zo {k) Dz E {k) 



(47) 
(48) 



Using the definitions of {Ik, Iq,Cq} and the proposed 
form of Es, Ea. l47|) transforms into 



-L=/ K ( 7 )(l-C Q ( 7) r)) 

'->n 

+I K (l + A 5 )(1 + C Q ( 7 + A s ,r)) , 
from which one obtains 

D 



7 (r) -De j ^- c q'--<)) ( 



A 5 + 7 



l + C Q (-, + A s ) 
) 1-Cq(t) 



(49) 



(50) 



where As = S — 5 e ; only the energy dependence of Cq is 
shown. This is the full expression for 7 when the odd dou- 
blet is the dominant one. For Cq(5k) = one has 5 — 
5 e = 8k, and Ea. (|oT7|l reduce to 7 2 = D 2 exp (— 1/J„), 
i.e. 7 = 5 K , see Ea. <|4T)|) . 

At the maxima of Cq (r ~ n\p) one has 5 3> <5 e ,7 , 
thus As ~ 5 , and using Eg. 1)35(1 one obtains 



7 = 1? exp {- 



1 



f 1 + C Q (6 ) 

Ml-Cgij)} 1 2 + C Q (5 ) 



\], (51) 



i.e. nearly Ea. ft-SOJ) . the previous "partial" result for 7 
at those regions. For these maxima, 7 is exponentially 
small if CQ("/,r) is still strong, as discussed in Section 

linn 

To analyze the r-regions in which Cg(5^,r) < one 
must propose Es — —2E c i — 5 e — 7 given that for these 
regions 5 < 5 e . The result for 7 is just Ea. H49|l with 
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FIG. 4: Correlation energy gains of the odd S and even 8 e 
Doublets and the Super-Singlet <5> + 7 as a function of the 
distance between the impurities, for J n — 0.07238 (5k = 
0.001D) in the ID case. The inset is 7, the additional energy 
gained in forming the Super-Singlet. For large r it tends to 
5k and the total energy of the Super-Singlet to —2Ed — 28k- 



the changes 8 <-» 5 e and Cq — > —Cq. In Fig.Q we 
plot the correlation energy gains of the Doublets and the 
Super-Singlet. 

At the first maxima of the doublets energy gains, i.e. 
the minima of their total energy, the additional energy 
gain of the Super-Singlet 7 is very small and therefore 
small thermal activation will destroy the correlation that 
generates it. This is not the case for the energy gain 
of the Doublets at those points, 5 a or 5 e , which is very 
strong, as discussed in Section UlI Bl 



C. The N T = 2 subspace 



In this subspace there are the lower energy configura- 
tions of the system, the AF singlet and the FM triplet, 
whose configurational energy is —2Ed- Given that these 
configurations are already the lower energy configura- 
tions of the subspace the effects of the hybridization on 
them can be directly evaluated by standard perturbation 
methods. For the sake of completeness we recalculated 
the result of reference 5 in the same variational wave func- 
tions framework as the previously analyzed subspaces. 
We use the "perturbation" property of the AF, FM states 
to solve their selfconsistent energy equation by series ex- 
pansion. For these AF and FM states two hybridization 
(H v ) steps are needed to bring out the RKKY energy 
term. For the S z — 1 component of the FM triplet, the 
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expanded VWF is 

|FM|) = |FM|) Z{k) 

k 

(5 fc 4 fcT 4 T + ^fe^ fcT 4 T )|F) + Y ( k > Q){ 



k , q 



(2A k A q 
+(A k A q b 

+i(S k A q bl k ^ Aq1 - A k S q b\ k ^ Sq] ) |AF) } . (52) 
After minimization, it turns out that 



y(fc, g ) = 

Z(k) = 



vZ(q) 
D Y (k,q) 
2v 



D z {k) 



(53) 
(54) 



9 \ 2 + cqr ckr , . 

E FM + E d -e k - 2v 2 n „. 3 ,(55) 



^z(fc) 

£fm = -2£ rf + 2v^Z(£;) , 



D Y (k,q) 



(56) 



where Dy{k, q) — E-pyi J r 2E d -'e k ~e q . As in the previous 
cases, the energy is to be determined selfconsistenly. We 
can make a Taylor series expansion on v 2 on the right of 
Eq.jnUJ), taking into account that Efm is also a function 
of v 2 . If this expansion is careful analyzed, it turns out 
that the expansion parameter is J n . The final result for 
-Efm, to order v , is the well known perturbation theory 
expression 



Evm — - —2Ed — v 



16 



^ E d + e k 

k 



k , q 



v> { 



16 + 8 cos g^r cos fc x r 
(£d + e fc ) 2 (£d + e,) ~ (£ d + e k ) 2 (e k + e q ) 



}, (57) 



the first v term in the rigth is the second order contribu- 
tion of the v 2 correction, this a well known property of 
variational wave functions like the Kondo Varma-Yafet 
singlet. In the second v 4 term there is a contribution 
that depends on the interimpurity distance r, 



v 4 ^ cos (q x r) cos (k x r) 

- 8 -pi 



= 4 In 2 D J 2 J R (r) 



(58) 



this is half the RKKY energy. The prefactor 8 appears 
because the sum in Ea. lj5Sll is over k x > (q x > 0) 
electron (hole) excitations. The r dependent part Jr{t) 
of the 1D-RKKY is given by 2 ^, 



Mr) 



-[- 

7T L 2 



Si(2fc F r)] 



(59) 



which is plotted in Fig.l. It decays very fast, its first ex- 
treme values are 1.0 at r = O.OAi?, and —0.179 at Af/4. 



At r = 2.0Ai? its value, 0.025, is already very small. The 
r = value of T, R is determined by integration of Eq. 
in the half filled flat band approximation for the con- 
duction band of the metallic host. The Kondo energy, 
Ecis. H5I6|) . is also evaluated in that approximation. For 
the AF state, expanded to the same order as the FM 
state, there appear the same terms as in Ea. H57f) but with 
a minus sign in front of X^, thus the RKKY interaction, 
i.e. the energy difference between the FM and AF states, 
is twice Efl. 



D. The RKKY and the Doublets 

The Doublets states are the results of a first-order in- 
teraction present in the TIA Hamiltonian. Only the con- 
figurations that are reached with one Hy step from the 
corresponding vertex state are included in the VWF we 
use to study them. See also that it depends on J n (al- 
though in a non-perturbative way) and that just one hole 
excitation is involved in the calculus of Cq . The RKKY 
is a second-order interaction, two Hy steps are needed in 
their VWF, it depends on J 2 and both a hole and an elec- 
tron excitation are involved in its evaluation. In fact, the 
correlation energy of the dominant Doublet is higher than 
the corresponding to the RKKY in most of the ( J n , r) pa- 
rameters space of the Hamiltonian. At r = the equa- 
tion £r(0) = S3 determines that for J n > 0.085 the 
Doublets correlation energy is greater than the RKKY 
energy. Note that the one-impurity Kondo energy 8k is 
never greater than X#(0). As a function of r the RKKY 
interaction decays at shorter distances than the Doublet 
one. In Fig0we show the "Doublets-RKKY" "Parame- 
ters Space Phase Diagram" for the ID two-impurity sys- 
tem, determined by max((5 ,<5 e ) = |Er|. This is not a 
true Quantum Phase diagram given that not the Dou- 
blets nor the FM (AF) states are the ground state of the 
system. But it indicates in what (J„,r) regions is rel- 
evant to include the second-order RKKY effects in the 
analysis of the Doublets structure. Over the near ver- 
tical dashed lines that mark the limits between the odd 
and even Doublets, determined by Cq(Sk) = 0, there is 
little enhancement of the Doublets energies. These zones 
coincide with the AF maxima of the RKKY interaction 
and thus over these lines there are the highest indents of 
the RKKY regions. 

How Efl(r) modifies the internal structure of the Dou- 
blets has been analyzed in Rcf. 24 , where the Doublets 
where first introduced. To include the RKKY effects 
in the VWF of the odd Doublet, two more Hy gen- 
erations of configurations must be added to it. In the 
odd Doublet VWF, Ea. (|31|l . see that the configurations 
with the ZA(k) amplitude factor correspond to a AF 
state of the impurities, whereas the ones with the Z$(k) 
factor are configurations with their impurity part in a 
FM state. As in Section ITVAl we cut our VWF at 
that order both amplitude factors result to be equal, 
Z A = Z s = v/(E Do - (-2E d + e k )) = -v/(6 a + e k ). 
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FIG. 5: Doublets-RKKY ID "Quantum Phase Diagram" for 
the TIA Hamiltonian. The horizontal dot lines mark the val- 
ues of &k (in units of D) for the corresponding values of J„. 



We show in Appendix IO that the main effect of the con- 
figurations needed to include the RKKY interaction in 
our Doublet VWF is to modify those amplitude factors. 
Their denominators will reflect now the different ener- 
gies of the FM and AF impurity configurations. The 
self-consistent equation for the odd Doublet energy when 
these higher order process are included is given by 



En = 



3v 2 (l + cos{k x r)) 



E D + 2E d + 2£ + \J X + S fl - e k 



v 2 (l — cos(k x r)) 



2Ed + 2£ — jJx — £fl — ek 



,(60) 



where the first sum in the right comes from the con- 
tribution of the fcrro configurations in the Doublet and 
the second one from the antiferro configurations, T. is 
the one impurity corrections (up to second order in J n ) 
and the terms proportional to Jx appears, just with the 
VWF used in Section IIV Al if an external RKKY-like 
term {—JxSl-Sr) is added to the TIA Hamiltonian. The 
effect of So is to shift the "zero" of energy and the ef- 
fect of Jx is similar to the one of the intrinsic RKKY 
term thus in what follows we discuss just the effects 
of this last term in the Doublets structure and energy. 
With the considerations above, the variational amplitude 
factor of the FM configurations in the Doublet is given 
by Z s (k) = -v/(A - S fl + e fc )) when RKKY effects 
are taken into account. And for the AF configurations 
Z A (k) = -v/(A + £ fl + e fe ) holds. A = |£ fl (r)| + <5 s (r) 
is the total correlation energy gain of the Doublet, com- 
posed by the RKKY and Kondo-Doublet contributions. 

For Y,r > (which corresponds also to the maxima of 
|Cq|), the dominant components of the dominant doublet 
are the FM configurations, thus we assume Ed = —2Ed — 
Y.r — 5f, and using the definitions of Ik and Cq, we 



(1 



1 

c Q (sf 



3(1 



com 

2 0? 



2Sh)) 



In 



D 



<E + 2S; 



(61) 



a self-consistent equation that easily gives the value of 
the odd-doublet correlation energy (5g) in a RKKY-FM 
region. In the extreme case Hr 3> S one can neglect 
the AF-configurations contributions to Ea. (|6ip (this is 
equivalent to use for the Doublets a VWF with Z A {k) = 
0) and thus 



D exp 



-1 



(3/2)(l + C Q )J n 



(62) 



which gives 5f > 8 k for values of \Cq\ down to 1/3. Note 
also that for high values of Cq there is little reduction of 
(5j compared with its bare (T,r = 0) value S , because the 
AF-configurations contribution (~ (1 — Cq)/2) is already 
very small in this case. 

Similarly, in the RKKY-AF regions (£r < 0, and 
Cq ~ 0), we assume Ed = —2Ed — |£r| — 5^ and, with 
the notation Ha = |£r| = —Hr, the general equation 
Eq.lJBUJl can be expressed as 



1 _ 3(1 + C Q (<5 S +2E^)) 



■J 



In- 



D 



(1 



-2S A 



(63) 



that in the extreme case XU 3> 5k can be reduced to 



5% = Dexp 



-1 



(i/2)(i-a Q )j n 



(64) 



which gives very small values for the correlation energy 
of the Doublet in this situation: for J n — 0.07238, which 
corresponds to 5 k = 0.001D and S 3 = S a (r = 0) = 0.01D, 
one obtains = 0.000000000001D. Therefore the for- 
mation of the Doublets in an extreme RKKY-AF region 
will be very difficult to detect in a numerical simulation of 
the system. Experimentally, nearly any thermal activity 
will destroy the doublets in this regime. 

The previous equations (Ecis. H60|) to l|64|l ) are for the 
odd Doublet, which dominates for Cq > 0. As the equa- 
tions corresponding to the even Doublet (which domi- 
nates for Cq < 0) are the same but with a minus in front 
of Cq, the final result is that Eas. (|61l64|l also apply to 
the even Doublet with the change Cq — > —Cq. In Fig^l 
we show a "generic" analysis of Ea. (|6l1|) : for J n = 0.07238 
and fixed values of Cq we plot 6^ as a function of Tir 
(— T<r) on the right (left) panel, energies are in units of 
8k- Actually, both Cq and Y>r are precise functions of 
(J„,r) and not any value of T,r can be obtained for a 
given J„. Anyway, we use this figure to show the general 
trend of the solutions of Ea. (|60|) . It can be seen that 
a positive ~Sr has little effect in the dominant Doublet 



11 




CO • ■ . , ■ , . — 4 4- — . , . , . 1-0 

-2.0 -1.5 -1.0 -0.5 0.0 5 10 15 

\ py 

FIG. 6: Correlation energy gain of the odd Doublet when a 
non-negligible RKKY interaction is present. Note the change 
of scales (and behavior) between the left panel, corresponding 
to the E_r < case, and the right panel, which corresponds 
to the Sjj > case. 



structure and it nearly has just an additive effect in the 
total energy of the Doublet. This response is due to the 
fact that the dominant Doublet is primary formed upon 
FM-like configurations. By the same reason a negative 
Eij, which energetically penalize the FM-like configura- 
tions, has a catastrophic effect on the energy gained by 
forming the Doublet. The corresponding extreme limits 
are well described by Eqs.|JB5J| and 
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contribution. Only at kpr = ir/2, a crossing point for 
the Doublets (Cq = 0), the RKKY-AF interaction pro- 
vides "all" the energy gain of the dominant Doublet. 
This point is an extreme RKKY-AF described in 

Eq.@. 

These extreme RKKY-AF regions, which for the actual 
TIA Hamiltonian are only present for low values of J n , 
are the subject of contradictory results in the literature. 
Numerical Renormalization Group (NRG) calculations 
seem to support a crossover from a Kondo- Singlet to an 
AF-Singlet for these regions. Instead, Quantum Monte 
Carloi^ and numerical VWF— calculations do not give 
support for such a quantum phase transition to occur in 
the U — > oo limit of the TIA Hamiltonian. Our analyt- 
ical VWF analysis is in line with these latter results for 
these extreme RKKY-AF cases: although with a very 
low correlation energy the Kondo-Doublets are formed 
on top of the AF configurations. Note that this result 
only applies in the analyzed limit of the TIA Hamilto- 
nian and assuming that the RKKY-AF effects can be 
analyzed in a perturbative way. We have use this last 
property to obtain the standard RKKY result in Sec- 
tion ^^21 and also in Appendix IU1 to obtain Ea. (|60[) . 
Note also that in this limit one has ^S> \Zs(k)\ 

( 1/(6? + e k ) » 1/(6? + 2|£ fl | + e k ) ) and thus the 
properties of the Doublet are determined mainly by its 
AF group of configurations. 

We have restricted our previous analysis to the sce- 
narios that are actually present in the analyzed limit of 
the TIA Hamiltonian: ferromagnetic RKKY regions co- 
incide with the maxima of |Cg(r)| and antifcrromagnetic 
RKKY regions are in correspondence with Cq(t) ~ 0. If 
one takes Cq and E/j as free parameters, as is often done 
in NRG simulations using the Kondo Hamiltonian, one 
can achieve the following situation: a strong AF RKKY 
and Cq = 1. In this case the self-consistent solution of 
Ea. l|64|) gives 6? = and the doublet is not formed. This 
is a possible explanation for the contradictory numerical 
results present in the literature. 



FIG. 7: Correlation energy of the Doublets including RKKY 
effects, for J n = 0.07238 and as a function of the inter- 
impurity distance r. The dominant Doublet energy is mainly 
composed by the Kondo-Doublet energy. Only at kpr = 7r/2 
the correlation energy is nearly fully provided by the RKKY 
interaction. This point corresponds to a maximum of the 
RKKY-AF energy and a zero of Cq . 

In Fig[7|we plot the total correlation energy of the Dou- 
blets (A x = 5 x + |E_r|, x = {o, e} ) as a function of the 
distance between the impurities and for J„ = 0.07238. As 
discussed previously the main contribution to the energy 
of the dominant Doublet comes from the Kondo-Doublet 
interaction, the RKKY contribution becomes negligible 
very rapidly when the distance between the impurities is 
increased beyond Af/2. For the dominated Doublet, in- 
stead, nearly all the energy gain comes from the RKKY 



V. Cq IN THE 2D AND 3D SCENARIOS 

In the previous Sections, when needed for graphical 
purposes, and because of its technological applications, 
we have evaluated our equations in the ID case. In ID 
the only decoherence factor, both for Cq and the RKKY, 
is the energy width of the packet of excitations involved 
in the interaction. In higher dimensions the angular de- 
coherence effect dominates the behavior of Cq at values 
of r lower than £k- In "sintetic" TIA systems (magnetic 
impurities deposited over a carbon-nanotube, Quantum 
Dots connected to a Quantum wire or to a confined 2D 
electron gas in a semiconductor heterostructure, etc.) 
the effective Dimension is generally within ID and 2D. 
The "classical" Two-Impurity Anderson problem corre- 
sponds to two magnetic impurities in a 3D metallic host. 
Here we analyze the behavior of 6 and 7 in the 2D and 
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3D Dimension cases. At this effect we only needed to 
evaluate Cq. For the 3D case it results 



r<3D 

and for the 2D case 



sin (kpr) /Uft , 



C^ D ~ BesselJ(0, k F r) , 



(65) 



(66) 



for r <C £,k, where the main decoherence effect is the 
angular one and they depend very weakly on 5. For r > 
£,k Cq decays like the power D (= 1, 2 or 3) of one over 
r. The RKKY Y,r has already the last behavior at r a 
fraction of Xf- 




FIG. 8: 2D Correlation energy gains 8 (dash line), 8 e (point 
line) and 8 > + 7 (full line) as a function of r, for values of J n 
corresponding to, from top to bottom, 8k = W~ 3 D, 10 -4 -D, 
and 10" 5 D. 

In Fig[H]we show three 2D cases corresponding, from 
top to bottom, to J„ = 0.07238, 0.05428 and 0.04343. 
The corresponding values of 8 k are 0.001 D, 0.0001 D, 
and 0.00001 D and the ones of 8 3 (—S o (0)), which not 
depend on the Dimension, are 10.0 8k, 21.5 8k, and 
46.4 8k, respectively. It can be seen that due to the 
angular decoherence effects the correlation energy gains 
decrease more quickly, as a function of r, than in the ID 
case (compare the top panel with FigHJ. It can also be 
seen that the relative effect of the Doublets interaction is 
higher the lower J„. 

In Fig[5]we show the 3D case for the same J„ values 
analyzed in Fig|Sj It can be seen the same trend than in 
going from ID to 2D. The top panel (3D, 8 K = 0.001D) is 
the case numerically analyzed in Rcf. [ 20 , Fig. 4(a)]. There 
are very minor differences, for example the numerical 
simulation does not reach the ^3 r — limit, although 
its authors already obtained that limit analytically. It is 




FIG. 9: 3D Correlation energy gains 8 (dash line), S e (point 
line) and <5> + 7 (full line) as a function of r, for the same 
three J n values analyzed in the 2D case. 



known that Kondo-like systems are very difficult to sim- 
ulate numerically due the logarithmic scales involved in 
the problem. 



VI. (Sl-Sr) correlation 

A very important observable in these systems is the 
impurities spin-spin correlation (Sl-Sr) because of its 
possible application in quantum computation. It is com- 
monly believed that this correlation is determined by the 
RKKY interaction. We have show in Section TlV Al that 
there is a simplest process than the RKKY one that leads 
to a impurity-impurity correlation, the Kondo-Doublets 
interaction. In Section HV Dl we have also shown that the 
correlation energy of the Doublets is greater than the 
RKKY correlation energy for most of the Hamiltonian 
parameter-space. Here we show that the Kondo-Doublets 
interaction generates a ferromagnetic alignment of the 
impurity spins without resort to the RKKY mechanism. 

Using for the Doublets VWF the simplest form given 
in Eq. l|31|l . i.e. without the inclusion of the higher or- 

the Doublet 



der RKKY effects discussed in Section llVDI 
impurity spin-spin correlation (Sl-Sr) is given by 



(D\S L -S R \D) 
(DID) 



= ± 



Dq(S, 



4 (2±D Q (8,r)) 



(67) 



where Dq — (8sIq) / '(OsIk) and the upper (lower) signs 
hold for the odd (even) Doublet. 

In FigEHit can be seen that the dominant Double fa- 
vor a parallel (ferromagnetic) alignment of the impurity 
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FIG. 10: Impurity-spin correlation for the Doublets, without 
higher order RKKY effects. The thicker section in each line 
is the r-region dominated by the corresponding Doublet. The 
straight segments at 1/4 (—3/4) corresponds to the RKKY 
prediction. 



spins. Note that this mainly ferromagnetic response is 
obtained without include the effects of the RKKY inter- 
action in the Doublets structure. 

This magnetic response can be clearly understand by a 
closer look to the first-order interaction behind the Dou- 
blets. Correlations between the impurities are the conse- 
quence of closed electron-paths that involve both impu- 
rities, paths driven by Hy in configurational space. For 
example, the second-order RKKY-FM path stars when 
one electron (say the one in the left impurity) hops to a 
band state k(> hp) and then a second electron hops into 
the impurity, letting a hole q(< kp) in the band. The 
path is closed by reversing these steps but at the other 
impurity. The matrix element for this loop, taking into 
account the mirror path, is 



m \h v \]]) R kky 



=Fv 4 cos (k + q).r 



(68) 



where the plus sign holds for the AF state. This correla- 
tion path produces the energy correction =pE/j(r) for the 
FM (AF) state, which splits the FM and AF states. We 
use in this Section the direct orbital representation used 
in the unsymmetrized Hamiltonian, Eq.Q). Therefore, 
in the equation above, |TT) = |FM|) = d L ^S R ^\F). 

In this representation the vertex state of the odd Dou- 
blet, d\^\F), is given by 

\A a ) ee i=(4 CT T d La )\F) ee -L(|0a) T |a0» , (69) 

where the plus sign is for the vertex of the even 
Doublet, the symmetric one-electron- in-the- two-impurity 
d Sa \F) = \S a ) state. The hybridization Hy connect the 
\A^) state, by promoting an electron from below kp to 
the empty impurity, with the following ones 

\A Uk ) = =1 4 T ( e —/ 2 |f!) +e +*-/2 | ir)) j(70a) 



V2 
— v 

71 



\A m ) = -J b{ L (e-^ 2 ITT) +e +Jk - r/2 |TT)) -(70b) 



Closing the loop, as depicted in Fig^] the following ma- 
trix elements are obtained 

(A 7 \H v \A Uk ) = v 2 , (71a) 
(A^HvlA^) = v 2 (l±cosk.r) , (71b) 

where the minus sign corresponds to a similar path but 
for the \Sf) state. The last element, which depends on 
the inter-impurity distance r, determines the properties 
of the Doublets. 

( " . t ) - ( T , - ) -E d ( - , T ) - ( T , - ) 

I t rxi 

( A , T ) ( T , 4 ) -2 E d + e k ( T , T ) ( T , t ) 



FIG. 11: First order loop paths for the d^lF) state. On the 
left it is shown the non-correlated aa channel. On the right 
it is the interference-enhanced aa channel, the crossed arrows 
correspond to the inter- impurity paths. 

As the starting state of the loop has a higher energy 
than the bunches of "visited" ones the energy associated 
with this interaction can not be obtained by perturbative 
methodsiiSS. This situation, a single state connected to 
bunches of lower energy states, is the hallmark of a Kondo 
structure, in this case the Kondo Doublets that we have 
analyzed in the previous Sections. 

With the analysis above at hand, the "ferromagnetic" 
behavior of the impurities induced by the Doublets is easy 
to understand. The correlated channel is the aa one, 
which corresponds to a ferromagnetic arrangement of 
the impurities (Eas. ij70bl71b|) ). The configurations con- 
nected via the non-correlated aa channel are also mainly 
ferromagnetic configurations at the maxima of 5d- See 
also Eas. (|33l I60|l . for the dominant Doublet the connec- 
tivity factor of the ferro configurations is ~ 3 ( 1 + 1 Cq | ) / 2 , 
whereas that of the antiferro ones is ~ (1 — |Cq|)/2. 
Clearly, the screening action of the hole is more effective 
when the impurity-spins are aligned and the distance be- 
tween them is near the resonant condition for the lowest 
energy holes (r ~ nXp/2, | cosfcpr| ~ 1), see Ea. (l71b|l . 

At the transition points from one Doublet to the other 
it is shown in FigE3a little negative value for (Sl-Sr) 
and then a jump to positive values. Actually, at those 
transition points the Super-Singlet has a maximum in its 
correlation energy (7 ~ 5 k) and it is formed with similar 
weights in both Doublets, thus an average of the response 
of the odd and even Doublets is to be expected at those 
regions, given a smooth transition with a (Sl-Sr) value 
near zero from one Doublet region to the next one. In 
FigEHit is shown also the RKKY prediction, based upon 
^fl( r ) ^ 0. Taking into account both interactions, the 
stronger of them determines the (Sl-Sr) response of the 
system, for max((5 ,<5 e ) ^ lE^I the response is deter- 
mined by the Doublets curve (Ea. l|o7|) . Fig ll(J|) whereas 
that in the opposite situation the impurity-spin correla- 
tion tends to that of the RKKY. 
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Now we include the RKKY effects in the Doublets wave 
functions, and thus in their (Sl.Sr) correlation. Follow- 
ing the analysis of Section IIVDI and Appendix El we 
take the variational amplitudes in the Doublets to be 
Zs ~ 1/(A_d — Tir + efc) for the ferro configurations and 
Za ~ l/(An + Si? + efc) for the antifcrro ones, where 
Ad = |Sfl| + <5|3 as evaluated in the RKKY-Doublets sec- 
tion. Thus, including RKKY effects, the impurity-spin 
correlation of the odd Doublet is given by 



(Sl-Sj- 



\ 3 Jf (1 



3 J A / 1 

4 AT I 1 



3 Jfc (1 



^(1 



where D 



S(A) 



J 



S(A) 

K 



D Q (A Do =f S fl ,r) and 



E 



n D A ( 



(72) 



(73) 



and for the even Doublet the change Dq i— > —Dq must 
be done. The supra S(A) terms come from the con- 
tribution of the ferro (antiferro) configurations of the 
Doublet. Eq.J72J) reduces to the previous Ea. (|u7Jl for 
5d{t) > |Sfl(r)|. Eq.{72J) reflects the changes in the 
energies of the FM and AF impurity configurations pro- 
duced by the RKKY interaction. Thus in a RKKY-FM 
region the relative variational amplitude of the ferro con- 
figurations is greater than the one of the antiferro ones 
+ e k) > 1/(2Er + Sd + efc), given a ferro-like 
(Sl-Sr) correlation. The opposite is true in a RKKY-AF 
region. Note instead that the ferromagnetic response in- 
duced by the Kondo-Doublet interaction depends on the 
interference enhanced matrix elements, the (1±Dq) fac- 
tors in Eq. l|72|l . not in the difference of energy between 
the FM and AF impurity configurations. 

In Fig^Jwe plot the < Sl-Sr > response of the Dou- 
blets, RKKY effects included, for the same ID case plotcd 
in Fig llOl There is an abrupt change in the response of 
the dominated Doublet (thin lines): as <5< is exponen- 
tially small (see Fig[3J) the RKKY effects determine the 
response of the dominated Doublet. Thus the thin lines 
closely follow the RKKY response, jumping from 1/4 to 
—3/4 (and back) at the £#(r) = points. Instead, the 
dominant Doublet response (thick sections) is little mod- 
ified by the RKKY interaction excepts at the extreme 
RKKY-AF point at k F r ~ 2, where \E R {r)\ > 5>(r) 
and thus the RKKY-AF behavior dominates. As dis- 
cussed in Fig0this point corresponds to a maximum of 
the antiferromagnetic RKKY and to a minimum of the 
Kondo-Doublet interaction (Cq ~ 0). Note also that the 
ferromagnetic response induced by the Kondo-Doublet 
interaction persists well into the RKKY-AF region. 

For the Super- Singlet, Ea. i(3*5jl . the (Sl-Sr) correlation 
is approximately given by a weighted average of the dom- 
inant and dominated Doublets response. The weighting 
factors are proportional to the average of the square of 
the corresponding amplitudes in the Super-Singlet VWF, 
Z [k) ~ l/(7 + e fc ) and Z E (k) ~ l/(j + A s + e k ) respec- 
tively (in the S > 8 e case), that turn out to be propor- 
tional to I/7 and 1/(7 + As). 
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FIG. 12: Impurity-spin correlation for the Doublets, RKKY 
effects included. The thicker section in each line is the r- 
region dominated by the corresponding Doublet. The domi- 
nated Doublet (thin lines) closely follows the RKKY predic- 
tion. 
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FIG. 13: (Sl.Sr) for the Super-Singlet (full line), the domi- 
nant Doublet (dash line) and the dominated Doublet (point 
line) as a function of the distance between the impurities. The 
upper panel is for a 3D case and the lower one for a 2D case. 



In Fig^H we plot the Super-Singlet and Doublets 
(Sl-Sr) correlation for a 3D and a 2D case (without 
RKKY effects) . See that both in the 2D and the 3D cases 
the ferromagnetic correlation induced by the Doublets is 
present for inter- impurity distances of several times . 



VII. CONCLUSIONS 

We have presented a simple and complete analysis of 
the two Anderson Impurity system. The variational wave 



15 



functions we use to analyze the problem are generated in 
each subspace of the Hamiltonian starting with the sim- 
plest configuration, and applying to it the hybridization 
terms, i.e. the non diagonal part of the Hamiltonian. The 
configurations so generated are grouped together with a 
variational amplitude function according with their in- 
ternal symmetry. This step is repeated as many times as 
needed to bring out the physics of the state. One step 
is necessary for the Doublet states, two for the Super- 
Singlet, one to connect the Doublets through the Fermi 
Sea state and one to retain the internal structure of the 
Doublets, two more are included in Appendix ICl to show 
how the Doublets are affected by the RKKY interaction. 
In this aspect, the technical method we use is in fact 
an analytical Lanczos method, to some point similar to 
the one used in Ref.f2£], but with variational functions 
as amplitude coefficients. Such variational functions are 
found by the classical method, i.e. Euler-Lagrange mini- 
mization of the energy functional, which is also evaluated 
analytically. 

This method can be also applied to the case of M im- 
purities, specially if the symmetry group of the impurity 
arrangement is contained in the symmetry group of the 
metallic host. In such a case M steps will be needed 
in the variational Af-supersinglet to arrive to the states 
with M electrons in the impurities, starting from the \F) 
vacuum state. But our conjecture is that, for a given 
set of system parameters the physical relevant state is 
the M-multiplet constructed by forming a Kondo singlet 
in the higher effective hybridization channel and filling 
the other effective impurities with one electron, as we do 
in Section IIII Bl to form our first approximation to the 
odd doublet. Starting from this M-multiplet, and going 
down, step after step, towards the |F) state, each succes- 
sive step will introduce a exponentially small correction 
to the M-multiplet energy, as 7 correct 8r> in the case we 
analyzed, Section IIII CI Other impurity configurations, 
like three in a row, are also of interest. The equations of 
the three-in-a-row case can be easily obtained from the 
present study, by adding an impurity at the origin. In 
Appendix [B] we outline this case and show that for cer- 
tain values of r the effective hybridization is very large 
in the active channel. A full analysis will be presented 
elsewhere. 

The TIA Kondo Doublet states generates a ferromag- 
netic Impurity-spin correlation, without resort to the 
RKKY interaction. The correlation energy gain of these 
Doublets, which is driven by a first-order direct co- 
herence effect of the hybridization, exceeds that of the 
second-order RKKY interaction for most of the Hamil- 
tonian parameter space (J„,r). These properties put 
the states we found well in the experimentally accessible 
range for Quantum Dots systems built on semiconduc- 
tors devices^, were most of the relevant parameters can 
be controlled by gate voltages. 

A milestone of the Two Magnetic Impurities problem is 
the renormalization group analysis of Ref. £»] . This work 
was carried out in the symmetrized basis too, but the r 



dependence of the effective hybridization was treated in 
a crude fashion, just taking its value at k = kjr, the rele- 
vance of the increased value of the effective hybridization 
was missed. Nevertheless, the two energy scales founded 
in that work can be traced to be the 7 and 5d correlation 
energy gains of the supersinglet and the doublets. Their 
coupling parameters J , J e are related to the square of 
our effective hybridization terms, i.e. ~ (1 ± Cq{5, r)), 
and their "Kondo temperatures" Tx(J e ), Tk{J ) to our 
Sd and 7, which are not the Kondo energies correspond- 
ing to J e and J but Ea.i|55jl and Ea. (|5T)|> . We show 
that their relative relevance depends on the value of 
Cq(Sk, r), for Cq{8k) — both are equally relevant and 
both tend to 5k- Our results are also in perfect accord 
with the VWF numerical analysis of L. C. Andreani and 
H. BeckSa. 

Summarizing, our results confirm the general picture 
of the system behavior already known from numerical 
simulations but, thanks to our explicit VWF that we an- 
alytically solve in the different subspaces of the TIA, they 
also point out important differences with the usual inter- 
pretation of that numerical data. The main interaction 
between the impurities, contrary to the general belief, is 
not generated by the RKKY process but by the interfer- 
ence enhanced hybridization that generates the Kondo- 
Doublet states. The strength of this Kondo-Doublet in- 
teraction depends on the one-hole quantum interference 
factor Cq{t). 

For inter-impurity distances r around the maxima of 
|CQ(r)| a two-stage Kondo screening takes place. It is 
a correlated quenching of the total spin of the system, 
not the successive "one-impurity" Kondo screening of the 
odd and even hybridization channel s 9 ! 15 . It also gener- 
ates a strong ferromagnetic correlation between the im- 
purity spins because the interference effects increase the 
weight of the ferro-like impurity configurations in the 
dominant Doublet. The ferromagnetic r-regions of the 
RKKY interaction are comprised in these zones. For 
low values of both the coupling constant J n and r the 
RKKY interaction is strongest than the Kondo-Doublet 
interaction (Er > <5> ^> 5k)- However this has little 
effect in the structure of the TIA Kondo states because 
the ferro-like impurity configurations are already their 
main component. In a "temperature scheme"— this sit- 
uation (Si? > <5>) generates the intermediate triplet spin 
1 phase between the high temperature uncorrelated im- 
purity spins phase and the Doublet spin 1/2 region. 

For r near the zeros of Cq(t) there is not enhance- 
ment of the Kondo-Doublet energies (5 ,5 e ~ 5k) and 
the system behavior tends to that of two uncorrelated 
Kondo impurities. The maxima of the antiferromagnetic 
RKKY coincide with these points. If — £r > 5k the 
antiferro-like impurity configurations dominate the struc- 
ture and response of the TIA Kondo states but, between 
the limits of our calculation, there is not a true quan- 
tum phase transition from the Kondo super-singlet to a 
AF-singletii. 

We provide detailed formulas for the evaluation, in all 
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dimensions and for any value of J„ and r, of the Kondo 
energies of the TIA Hamiltonian. We provide also ex- 
plicit formulas for the evaluation of the impurity spin- 
spin correlation in a given experimental situation. As we 
provide the explicit VWF of the involved states, other 
correlations can also be easily evaluated. 



origin and the third coefficient is the one correspond- 
ing to the Right impurity. The hybridization matrix el- 
ements, with the symmetric (%/2 cos fc^x) and antisym- 
metric (i\/2sinkxx) conduction band orbitals are 

V A = (d A \H hyb \c Ak ) = 2iv smk x r/2 , (B4) 
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APPENDIX A: M-CHANNEL SINGLET 

Given that for an arrangement of M magnetic im- 
purities is always possible to construct an M -channel 
singlet, the first step of the M-supersinglet, is worth- 
while to show that, although 8%m 3* $k, its total energy, 
Em = —Ed — 52M is always greater than the energy of M 
decoupled Kondo singlets, Ed c = M(—Ed — 8k)- There- 
fore, Eq.J^l), we must show that 



J2M 



-2M v 2 



1 



$2M 



(Al) 



is always greater than 



-M(Ed + 5k) = —M 2 v 2 - — ^ . (A2) 



Comparing the sums term to term, we see that in each 



for So M > 5 K . Thus E M > E, 



dc 



of them -j~ — . — > t — - 

OK+e q 6lM 

independently of the D, Ed 3> 5 usual approximation and 
for any distribution of the conduction band hole states 
(q). In particular, this applies for M= 2, Section UlI Al 



APPENDIX B: THREE-IN-A-ROW 

For the three-in-a-row impurity case, two as analyzed 
in this paper plus one at the origin, one can take as nor- 
mal modes of the impurities the following symmetrized 
combinations of the impurity orbitals 

A=(-l/V2, , 1/V2) , (Bl) 

51 - (1/V3, l/\/3, l/\/3) , (B2) 

5 2 = (1/V6, -2/V6, l/Vo) , (B3) 

In each of these combinations the first coefficient is the 
amplitude of the normal mode at the Left impurity, the 
second coefficient is the amplitude at the impurity at the 



V Sl = (d S i\H hyb \c S k) = \l - v (1 + 2cosk x r/2) ,(B5) 



Vs 2 = (d S 2\H hvb \csk) = \/g v (-2 + 2 cos k x r/ '2) (B6) 



for the symmetric combinations S. Therefore for r = 2\p 
one has {Va — 0, Vs 1 ~ V6 v, Vs 2 — 0} and a strong 3- 
multiplet can be proposed, formed upon a Kondo singlet 
in the Si impurity plus one electron in each one of the 
other impurity normal modes. 



APPENDIX C: THE ODD DOUBLET IN AN 
EXTREME RKKY-AF REGION 

In order to show the validity of Eq. (|60|l to analyze the 
RKKY effects in our doublets, we present here a partial 
deduction of it. Assuming we are in an extreme RKKY- 
AF region, we can drop the contribution of the FM-like 
impurity configurations to the odd doublet. Therefore 
we start with 



D 



of, 



= d^\F)+iJ2z A (k) A k 6 f 



Ak] 



|AF), (CI) 



applying Hy to the Z A configurations one electron is re- 
moved from the impurities and transferred to the band 
states. Given the structure of |AF) there are four possi- 
bilities, two corresponding to the transfer of a symmetric 
electron (with a probability amplitude proportional to 
S q ) and the other two corresponding to the transfer of 
an antisymmetric electron (~ iA q ), therefore the follow- 
ing states are added to the VWF 



E A k 



Ak] 



{ 



k.q 



i S q [Yi(fc,<7) 4 qL dl r \F)-Y 2 (k,q) 4^4jF)] - 
A q [Y 3 (k,q) c AqL d A1 \F)-Y 4 (k,q) c\^d\ { |F)]},(C2) 

applying again Hy, a big set of new configurations is 
generated. The ones with one electron in each impurity 
are the ones of lower energy, and of these, in order to 
maintain this presentation as simple as possible, we in- 
clude here only the ones that contribute to the RKKY 
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interaction 



k,q,p 



S q A p [Xi(k,q,p) b 



t J 

Apl C Sql 



|FMt) - 
X 2 (k,q,p) ^ pT 4 T |FM|>] + 
[X 3 (k,q,p) bl pl c\ ql |FMT>- 
X 4 (fc,g,p)6t pTC ^ T |FMi)]}. (C3) 

The discarded ones, with their impurity part in a |FM0) 
or a |AF) state, contribute only to the "one impurity" 
correction. They were included in Section ITvTI The 
energy of the Doublet is given by 



E D = (D oT |ff|D oT )/(D oT |D oT 



(C4) 



The variation of Eq. l|C4|) with respect of Z A (k), , 

Xi{k,q,p) gives nine coupled equation. These equations 
can be used to reduce Ea. ilC4(l to the form 



(C5) 



The variational equations are solved in a progressive way, 
starting with the higher order amplitude factors 

Xx{k,q,p) = -2vY 1 (k,q)/D x (k,q,p), 

X 2 (k,q,p) = -v (Y 2 (k,q)-Y 2 (p,q))/D x (k,q,p), 

X 3 (k,q,p) = -2vY 3 (k,q)/D x (k,q,p), 

X A (k,q,p) = -2vY 4 (k,q)/D x (k,q,p), (C6) 

where 



D x (k, q,p) = -E D - 2E d + e k + e q + e p . (C7) 

Using the above results, the Yi factors are found 

Y 1 (k,q) = -2v Z A (k)/D YS (k,q) , 
Y 2 (k,q) = -2v Z A (k)/D YS (k,q) + 

4v 2 ^(A 2 Y 2 (p, q)/D x (k, q,p))D YS {k, q) , 
p 

Y 3 {k,q) = -2v Z A {k)/D YC (k,q) , 

Y 4 (k,q) = -2v Z A (k)/D YC (k,q) , (C8) 



where 



D YS (k, q) = -E D - E d + e k + e q 

-4v 2 J2A 2 p /D x (k,q,p) , 
p 

D Y c(k, q) = -E D - E d + e k + e q 
-4v 2 J2Sp/D x (k,q,p) . 



(C9) 



The Yi factors carry the RKKY effects into the varia- 
tional equation obtained for Z A , which is 



(E D + 2E d - e k )Z A {k) = v 
2 



-v^{cos 2 (^)(y 1 (fc, g ) + y 2 (fc, (? )) 



+ S m 2 (^)(Y 3 (k,q)+Y 4 (k,q))} . (CIO) 

This is a self-consistent equation, the Yi factors depend 
themselves on Z A . In order to solve it, and to obtain the 
standard expression for the RKKY contribution, the Yi 
factors must be expanded in powers of v, as already done 
in Section IrvTl 

Yi(k,q) = Z A (k)[vy (k,q)+v 3 y 2S (k,q)} , 

Y 2 (k,q) = Y 1 (k,q)+v 3 y 2X (k,q) , 

Y 3 (k,q) = Z A (k)[wy (k iq )+v 3 y 2C (k,q)l 

Y 4 (k,q) = Y 3 (k,q) , (Gil) 

where, with D Y (k, q) = —Ed — E d + e k + e q , 



y (k,q) = -2/D Y (k,q) 

-8 \ - . 2 f px_r 
2 



V2s(k,q) = 



D Y (k,q) 



J^sin 2 (-^-)/D x (k,q,p), 



V2c{k,q) 
V2x(k,q) 



-8 ^ Z A (p) sin 2 ( ) 
D Y (k,q) ^ D Y (q,p)D x (k,q,p)- ( ' 



Using the above equations in the variational equation for 
Z A (k), it is obtained 



Z A (k) 



-v(i - v 3 r; 



-E D - 2E d + e k - S z ' 



(C13) 



where 



Tz = ^ ^( P )sin 2 (g|I)cos 2 (2|I) 



D Y {k,q)D Y (q,p)D x (k,q,py 



YZ Z = 4v 2 V 

4- D Y (k,q) 



8v 4 V- 
^ D 



q.p 



1 - cos(q x r) cos(p x r) 
(k,q)D x (k,q,p) 



2£ i + S 2 - S fl (r) . 



(C14) 



Both Tz and Y*z are very weakly dependent on e k , in 
fact it is a common practice to substitute factors of the 
kind l/{E d + e k ) by 1/E d in the available calculus of the 
RKKY interaction, a notable exception is ref.fSZ] . In 
any case, such approximation overestimates the RKKY 
effects. Soi is the order J % n "one impurity" correction. Tz 
can be evaluated by substituting Z A (p) by its first-order 
approximation v/ (Ed + 2E d — e p ), Tz gives a J 2 correc- 
tion to the "connectivity" factor of the AF configurations 
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involved in the odd doublet. In the equations above the 
sums over k,p (g) are over "symmetrized" hole (electron) 
excitations. From Eas. l|C5IC13l and IC14| the last term 
in the right side of Ea. (|6U|l immediately follows. 

We have analyzed the general situation, i.e. includ- 
ing the Zs FM-like configurations in the doublet (and 
the corresponding secondary configurations), in order 



to fully certify Ea. (|60() . Our "second quantization" 
Mathcmatica® package does the calculations needed for 
Eq. (|C4|I on the fly, a primitive version of this package was 
used in Refs.fS&Si]. To bring those equations to ink, in- 
stead, will take a considerable amount of pages and time. 
The general procedure is similar to the one we outlined 
here for the extreme RKKY-AF case. 
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